Scaler mode of the Auger Observatory and Sunspots 



Carlos A. Garcia Canal 
Instituto de Fisica La Plata, CCT La Plata, CONICET and 
Departamento de Fisica, Facultad de Ciencias Exactas, 
Universidad Nacional de La Plata 

a; 

^ ■ CC 67, 1900 La Plata, Argentina 

Carlos Ho j vat 

Fermilab, P.O. Box 500 - Batavia IL 60510-0500, USA. 

Tatiana Tarutina 
Instituto de Fisica La Plata, CCT La Plata, CONICET and 
Departamento de Fisica, Facultad de Ciencias Exactas, 
Universidad Nacional de La Plata 
CC 67, 1900 La Plata, Argentina 

Received ; accepted 



- 2 - 



Abstract 



Recent data from the Auger Observatory on low energy secondary cosmic ray 
particles are analyzed to study temporal correlations together with data on the 
daily sunspot numbers and neutron monitor data. Standard spectral analysis 
demonstrates that the available data shows 1/ f 13 fluctuations with f3 1 in the 
low frequency range. All data behave like Brownian fluctuations in the high fre- 
quency range. The existence of long-range correlations in the data was confirmed 
by detrended fluctuation analysis. The real data confirmed the correlation be- 
tween the scaling exponent of the detrended analysis and the exponent of the 
spectral analysis. 

Subject headings: Auger scalers, sunspots, 1/f fluctuations, scaling exponent 
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Introduction 



Solar activity gives rise to a modulation of the flux o: 



cosmic rays observed at Earth. 



The Pierre Auger Observatory (IPierre Auger Observatory! ) has made available the scaler 
singles rates observed on their surface detectors reflecting the counting rates of low energ y 



secondary cosmic ray particles ( IPierre Auger Collaboration 



2008 



2011 



Dasso et al. 



20121 ) . 



These data are presented after corrections for atmospheric effects, pressure in particular, 
and compared wit h temporal variations of solar act ivity as shown by data obtained with 



neutron monitors ( INeutron Monitor Database 



201 ll ). Solar and cosmic rays data can be 



presented as a temporal series containing modulations, correlations and noise fluctuations. 



'he availabi 



(ISIDC-team 



ityof Auger scaler data and data from neutron monitors and sunspot numbers 



201 ll ) motivated us to study the existence of long-range correlations present in 



the corresponding time series. 

The Pierre Auger Observatory is located in the city of Malarge, Mendoza, Argentina. In 
order to study the highest energy cosmic rays it covers a surface of 3,000 square kilometers 
with 1600 surface detectors sensitive to the transversal distributions of the showers 
generated by the primary cosmic rays. In addition, the longitudinal shower development is 
measured by atmospheric fluorescence by 27 optical telescopes. To monitor performance, 
the surface detectors have scalers that record signals received above a given threshold 
independent of any further shower reconstruction involving other neighboring detectors. 
This is refer to as the "scaler mode" of the Auger surface detector. The low threshold 
rates or scaler data have been recorded by the surface detectors of Auger Observatory since 
March 2005. These data should be sensitive to transient events such as Gamma Ray Bursts 
and solar flares. The rates at each detecto r are registered every second a nd the 15 minute 



average rates are available for public use ( lAuger scaler data online 



20111 ) . The temporal 



variations can be accurately studied, as these rates are very large as compared to other 
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data on solar activity. 



In the work of iPierre Auger Collaboration! (120111 ) t he pressure co rrected Auger scalers 



were compared to data from the Rome neutron monitor (IStorini et al. 



) and it was concluded 



that Auger scalers could be suitable for the study of solar activity. 

A sunspot is a temporary phenomenon in the solar photosphere t hat appears like a 



dark y isible spot compared to the surrounding regions (see, for example, 



Bray &: Loughhead 



( 119791 )). It corresponds to a relatively cool area of the Sun photosphere (1500 K less than 
the average photosphere temperature) as a result of the heat convection process inhibition 
by intense magnetic fields. 

The number of sunspots and their position on the Sun face change with time as a 
consequence of the solar activity cycle. The maximum solar activity corresponds to a large 
number of sunspots and in the minimum less sunspots are observed. The spots usually 
appear in groups. 

Data on sunspots are available for the last four centuries. From 1749 to 1981 the 
sunspot data were provided by the Zurich Observatory. Nowadays the World Data Center 
for the Sunspot Index in the Royal Obser vatory of Belgium is responsible for recording 



sunspot data. All data is available online ( jSIDC-team 



20111 ). 



Neutron monitors are ground based detectors that measure the flux of cosmic rays from 
the Sun and low-energy cosmic rays from elsewhere in the Universe. In a typical neutron 
monitor, low-energy neutrons produced by nuclear reactions in lead are slowed down to 
thermal energies by a moderator and detected by proportional counter tubes. A worldwide 
network consisting of approximately 50 s tations is in operation and their data are available 



on-line ( jNeutron Monitor Database 



201lh . 



For this analysis we select data from two neutron monitoring stations. Those sites were 
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chosen by the availability of complete data for all the days of the period of availability of 
Auger scaler data. Therefore there was no nece ssity to perform an int erpolati on procedure. 



These 2 neutron monitor are known as JUNG ((jJUNGj)) and APTY ( fjAPTYl )). The JUNG 
detector is located on top of the Sphinx Observatory Jungfraujoch, Switzerland and APTY 
is situated in the town of Apatity, Russia. 

The aim of this note is to present a systematic analysis of the temporal series from 
different experimental determinations. This analysis, b ased on power spectra behavior 



( Fan chiotti et al. 



1994 ; iFanchiotti et al. 



20041 ) and detrended power behavior ( IPeng et al 
20041 ) allows us to gain quantitative information about correlations among the different 
phenomena. It is also interesting to study the connection between the information provided 
by the power spectra analysis and the detrended analysis for the case of real data. 

This paper is organized as follows: in section 2 we present the results of the power 
spectra analysis of the three data sets and in section 3 we present the results of the 
corresponding detrended analysis. In section 4 final remarks are given. 



2. Power Spectra 



There is a number of accepted techniques ([Hamilton! 1 19941 ) to characterize random 
processes x(t). Usually one starts with the correlation function defined as 



G(t) = (x{t )x(t + r)) t0 - (x(t )) 



Another widely used tool is the frequency spectrum that is defined as the 
squared amplitude of the Fourier transform of the time signal: 



S(f) = lim - 

T— ¥oo 1 



drx(t)e 



2infr 



(2) 
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For a stationary process, the frequency spectrum is con nected with the 



correlation function through the Wiener- Khint chine relation (IMacDonald 



;emporal 



19621 ): 



S(f) = 2 drG(r) cos(2n fr). 
Jo 



(3) 



A true random process, also called white noise has no correlations in time, therefore 
the correlation function for the white noise is a delta-fu nction and the spectrum function 



19621 ) or random walk corresponds to the 



S(f) oc const. The Brownian motion (IMacDonald 
spectrum function S(f) oc l// 2 . 

Many naturally occurring fluctuations of physical, astronomical, biological, economic, 



traffic and musica l quantities exhibit S ( f) oc 



(see, fo r 



Bailliel jl996j ); 



example. IVoss fc Clar 



Novikov et al. 



ke! dl97fih : 



Press 



f be h avior over a l 



( 119781 ): Matthaeusl ( 119861 ): 



measured time scales 



Van Vliet 



(1199lh : 



( 1l997l )). These fluctuations are of interest because they 



correspond to the exist ence of e xtrem ely long-range time correlations in the time signals. 



This can be shown (see iJensenl (119981 )) if one assumes that the spectral function of a time 
signal S(f) ccl/f 13 and that the temporal correlation function G(t) oc l/r a . It follows from 
([3]) that 1/f 13 ex l/r 1_a . When m 1 it follows that a ~ which corresponds to correlation 
function G(t) ~ 1. 



An analysis of monthly sunspot data was performed in the work of iFanchiotti et al. 



(120041 ) . including the study of the power spectra and detrended analysis. It was found that 
the high frequency part of the spectral function of the monthly sunspots shows the 1/f 13 
behavior with = 0.8 ± 0.2. This corresponds to the 1/f noise. 

In Figure [1] we present the available ti me series data obtained as in dicated in 



corresponding refe rences: (1) Auger scalers (jAuger scaler data on line 



(ISIDC-team 



20111 ) and (3) JUNG neutron monitor data (jNeutron Monitor Database 



201lh. (2) sunspots 



2011 



) 
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for the time period when data on the Auger scalers are available. As the Auger scalers 
contained gaps, sometimes of various days, in order to use it in the analysis we performed an 
interpolation of the data justified by the Brownian behavior of data for the high frequencies. 

In Figure [2] we present the power spectra of available data on a logarithmic scale vs. 
frequency. The power spectra for the sunspots and the neutron monitor were shifted with 
respect to each other to avoid overlap. The spectral function indicates the presence of two 
different behaviors for low and high frequencies that can be approximated by linear function 
with the smaller slope in the range of lower frequencies. The frequency corresponding to 
the change from one region to another is different for each data set. We also note that in 
the sunspot daily spectra there is a peak corresponding to the frequency of approximately 
1/27=0.04 days -1 , not clearly seen on Figure [2] because of the low statistics. This 27-day 



peak is due to solar rotation (jBray fc Loughhead 



19791 ) . At first sight the description 



with the piecewise linear function seems more justified for the spectral function of the 
sunspots than in the spectral function of the scalers which shows more attenuated behavior 
for relatively small frequencies. In addition, the slope for the low frequency part of the 
spectrum for the Auger scalers and the neutron monitor is smaller compared with that of 
the sunspots. We make a further analysis of these facts below. 

We analyzed the power spectra of Auger scalers, sunspots and two neutron monitors 
using a linear least squares fit. First we fitted the data for the total range of frequencies. 
The resulting /3 is presented in the third line of Table 1. All data have similar slopes when 
fitted over the total range of frequencies giving /3 ~ 1.6. 

Second, we performed a fit of the data to the piecewise linear function consisting of 
two parts with different slopes. The frequency corresponding to the point of change of the 
slope f c was adjusted independently for each data set to obtain the best fit and is presented 
in the second line of Table 1. The fitted values of /3 for the range of low frequencies are 
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daily average Auger scaler rate 




2006 2007 2008 2009 2010 2011 

daily average NM JUNG rate 




Date 



Fig. 1. — (Color online) Auger scalers, sunspots and neutron monitor JUNG data. 




10~ 3 10~ 2 10~ 1 

f [1/day] 



Fig. 2. — (Color online) Power spectra of daily Auger scalers, sunspots, and the neutron 
monitor JUNG. 
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• sunspots 



icr 3 icr 2 icr 1 

f [1/day] 

Fig. 3. — (Color online) Power spectra of daily Auger scalers and sunspots with the corre- 
sponding result of the fit with piecewise linear function (solid lines). Dashed lines indicate 
95% prediction bounds. 
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given in the forth line and for high frequencies - in the fifth line. In Figure [3] we present the 
results of the fitting for Auger scalers and sunspots. 

It is seen from this analysis that in the low frequency range all data show the existence 
of long-range correlations with « 1. For high frequencies, the spectral function shows the 
Brownian behavior. It is also seen that the position of f c is different for each the data set. 
It should be noted that in the case of sunspot s the position of f r could be dictated by the 



existence of 27-day cycle of the Solar activity (IBray fc Loughhead 



1979|). 



(a) 
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f. [1/day] 



Fig. 4. — The dependence of the power spectra slope (3 on the range of the frequency used 
for the fit: (a) low frequency range; (b) high frequency range. See text for explanation. The 
dashed lines are to guide the eye. 

In order to have more quantitative conclusions about the value of f3 we performed a 
systematic fit of all data for different frequency ranges. In this analysis we set various sizes 
of the considered frequency intervals [10 -3 day -1 , /J with = 10~ 3 day -1 + Af and up to 
approximately f c with Af = 10 0,2 day -1 and present the fitted slope (3 for sunspots, Auger 
scalers and neutron monitor JUNG in Figure H^a). For f] w 10 _1 day~ x , that is when fitting 
all sets of data in the range of low frequencies, namely 10~ 3 day" 1 < < f c , (3 ~ 1. In the 
case of sunspots for / m 0.04 day -1 the drop in (3 can be explained by the influence of the 
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mentioned 27-day peak which is not seen in the case of Auger scalers and neutron monitors. 
For smaller frequencies the error is large for all the data because of the low statistics but 
the results agree within the error bars. 

In Figure 11(b) we show the results of a similar analysis but this time for the range 
of high frequencies. Here we obtain the value of /3 fitting the spectral function over the 
frequency interval with varying size [fj, Jn], where fj = Jn — Af with Af as before, and 
J'n stands for the Nyquist frequency. It is seen that for fj « lCT 1 day -1 all the data predict 
a similar value of 0. As fj approaches 10~ 2 day -1 its value drops because of the influence 
of the part of the low frequency spectra. 

For the high frequency range, the data show agreement with the following behavior: 

S(f) « /- 19±0 - 2 (4) 

which corresponds to Brownian fluctuations. This is an expected behavior because at high 
frequencies the fluctuations become more random or less correlated. 

We can conclude that the data show the coexistence of two behaviors in the power 
spectra: (1) consistent with a 1/ f dependence for low frequencies (2) consistent with a 1/ f 2 



behavior for high frequencies. Our analysis confirms the result of 



Fanchiotti et all (120041 ) 



where monthly data on the sunspots was analyzed. The su nspots in this frequenc y region 



have 1/f behavior. The 1/f 2 behavior was not detected in iFanchiotti et al.l ( 120041 ) because 



it corresponds to the high frequencies which are not present in monthly data. 



3. Detrended fluctuation analysis 



Another statistical method to reveal the extent of long-range temporal correlations in 
time series i s the Detrended Fluctuation Analysis (DFA). This method was introduced in 
Peng et al.l (Il994[ ) and applied in many areas of research, including physical and biological 
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sciences. 



Consider the time series x(t) consisting of N samples. The procedure is as follows: 
1. First generate a new time series u(t n ) by means of: 

n 

i(t n ) = ^2x(t k ), u(0) = 0, n = l,...,N. 



a 



(5) 



k=i 



2. Divide the time series u(t) into non-overlapping intervals of equal length s. In 
each interval the data is fitted using least-squares to t he first order polynomial (the 



polynomial of 2nd, 3d order may be used as well, see 
giving a local trend y(t) = at + b. 



Kantelhardt et al. 



(j2002|)) 



3. In each interval calculate the detrended fluctuation function using: 

,N 



(k+l)s 

[F(s)] 2 = [*(*)- v(*)] 2 ; k = o,i 

t=ks+l 



4. Calculate the average of F(s) over N/s intervals. 

For a fluctuating time series, the expected behavior is as follows: 

(F(s)) ~ s Ha 



quantity 


Auger scalers 


Sunspots 


NM JUNG 


NM APTY 


log^/day- 1 ) 


-0.87 


-1.35 


-1.25 


-1.15 


(3 total / range 
P low / 
/3 high / 


1.525 ± 0.073 
1.28 ± 0.14 
1.91 ± 0.23 


1.591 ± 0.080 
0.55 ± 0.32 
2.02 ± 0.13 


1.564 ± 0.078 
1.12 ± 0.23 
1.79 ± 0.14 


1.549 ± 0.079 
1.11 ± 0.22 
1.86 ± 0.16 



Table 1: Values of f3 for different frequency ranges obtained in fitting the power spectra with 
piecewise linear function with f c as a free parameter in the fit. 



13 



wher e H Q is a scaling exponent, a generalization of the Hurst exponent ( IKantelhardt et al 



20021 ). 



Initially the detrended fluctuation analysis was proposed as an independent measure of 
long-term correl ation, complementary to the spectral analysis information. In the work of 



Buldyrevl (119951 ) it was noted that the value of ft of the spectral function and ft' = 2H a — 1 
are "remarkably close to each other". In this analysis H a = 1/2 corresponds to the case of 
the white noise, 1/2 < H a < 1 reveals the existence of positive correlations and the special 
case when H n = 1 correspond s to 1/f noise. Random walk is characterized by H a = 3/2. 



Heneghan Sz McDarbyl (120001 ) examined the analytical link between DFA and spectral 
analysis and showed that they are related through an integral transform. It was concluded 
that DFA and spectral measures provide equivalent characterizations of stochastic signals 
with long-term correlation. In this work we study the connection between the scaling 
exponent H a and the coefficient ft introduced in the spectral function analysis to check this 
assertion in this particular case with the real data. 

In Figure the result of the detrended fluctuation analysis of Auger scalers, sunspot 
numbers and one neutron monitor (JUNG) is shown. In the case of the sunspots, the 
data can be separated into 3 regions with different slopes: (a) s < s c ; (b) s c < s < 10 2 
and s > 10 2 (the numbers here are approximate). The value of s c in the case of sunspots 
approximately corresponds to the 27-day peak in the spectral function. For the case of 
Auger scalers and neutron monitors the existence of three regions in the DFA curve is not 
as clear as in the case of the sunspots but a close examination of Figure \5\ suggests a change 
of slope at s ~ 10 for Auger scalers and neutron monitors. The situation is similar to the 
power spectra curves where the difference between the regions with 1/f and 1/f 2 is more 
apparent for the case of the sunspots. There is a correspondence between the frequency 
in the power spectra of the time series and the value of s. The region of large values of s 
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corresponds to the very small frequencies of the time series and have large fluctuations that 
can be seen on the Figure This region of large values of s is excluded from our analysis. 

In Table [2] we compare the exponent (3, and with the coefficient H a which 
characterized the slope obtained in the DFA analysis for two ranges of the frequencies: 
(a) "high frequency" s < s c and (b) "low frequency" s c < s < 10 2 , where log(s c ) is taken 
equal to the modulus of the logarithm of the frequency corresponding to the change of 
the regimen in the power spectra (and given in Tabled]). First of all it is seen that the 
calculated scaling exponent for low frequencies is always 1/2 < H a < 1 which implies 
the existence of long-range correlations. It should be noted that the analysis of the daily 
sunspot num bers confirmed the resu lt of the analysis of the monthly sunspot number 



performed by 



Fanchiotti et al. 



(120041 ). where the value obtained for the scaling exponent 
was H a = 0.62 ± 0.4. It corresponds to the low frequencies as the high frequencies are not 
seen in the monthly data. 

In order to compare the results of power spectra analysis and DFA we present in 
Figure M the slopes calculated using these two methods with the error bars indicated. The 
results are grouped by frequency: to the left for low frequencies and to the right for high 
frequencies. In each of these groups we present the resulting 2±1 from spectral analysis 
(filled markers) and the scaling exponent from DFA (open markers). It is seen that the 
results of both methods agree within the error bars and are consistent with H a w 1 for 
low frequencies and H a « 3/2 for high frequencies. Thus, in this work using real data of 
three different sources we confirm the correspondence between the value of /3 of the spectral 
function and the scaling exponent of detrended fluctuation analysis for both frequency 
ranges. 
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frequency 




0+1 
2 


H a 


Sunspots 








high 


2.02 ±0.13 


1.51 ±0.08 


1.40 ±0.06 


low 


0.55 ±0.32 


0.78 ±0.16 


0.75 ±0.03 


Auger scalers 








high 


1.91 ±0.23 


1.46 ±0.12 


1.39 ±0.16 


low 


1.28 ±0.14 


1.14 ±0.07 


1.14 ±0.03 


NM JUNG 








high 


1.79 ±0.14 


1.40 ±0.07 


1.33 ±0.09 


low 


1.12 ±0.23 


1.06 ±0.12 


0.96 ±0.04 


NM APTY 








high 


1.86 ±0.16 


1.43 ±0.08 


1.38 ±0.11 


low 


1.11 ±0.22 


1.06 ±0.06 


0.97 ±0.03 



Table 2: The result of the spectra analysis &±- and the scaling exponent H a from DFA for 
sunspots, Auger scalers and two neutron monitors. 
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4. Final remarks 

We have presented a statistical study of the available data on Auger scalers, sunspots 
and neutron monitors. We studied the frequency spectra and performed a detrended 
analysis. We found that the spectral function can be separated into two regions: (1) the 
region of the low frequencies and (2) the region of high frequencies that can be described by 
power laws. It was shown that the low-frequency part of the spectral function of all data 
shows 1/f behavior. Because of the low statistics the error is large. The high frequency 
part of the spectral function of all data is shown to behave like Brownian fluctuations. This 
similar behavior of the time series analyzed can be understood because all three kinds of 
events are correlated since they have their origin in the solar activity. 

The detrended analysis performed confirmed the existence of long-range correlations, 
for low frequencies, in all available data (the scaling coefficient H a 1). Finally, the 
correspondence between the scaling exponent and the j3 exponent of the spectral analysis 
was confirmed with real data. 
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Fig. 5. — (Color online) Detrended fluctuation analysis curves of Auger scalers, sunspots and 
neutron monitor JUNG. 
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Fig. 6. — (Color online) Comparison between and the scaling exponent for sunspot, 
Auger scaler and neutron monitor data calculated in the regions of high and low frequencies. 



- 19 - 



REFERENCES 



Auger collaboration 2004, Nucl. Ins. Meth. 2004, A 523, 50 Pierre Auger Observatory, 



http: / / www.auger.org 



Pierre Auger Collaboration, Bertou X. 2008, in Proceedings of the 30th International 
Cosmic Ray Conference (ICRC 2007) (Mexico DF, Mexico), pg.441-444 

Pierre Auger Observatory Collaboration 2011, JINST, 6, P01003 

Dasso, S., Asorey, H., For The Pierre Auger Collaboration, Advances in Space Research 
2012, 49, 1563 



Neutron Monitor Database, http://www.nmdb.eu/nest/search.php 



SIDC-team, World Data Center for the Sunspot Index, Royal Observatory of Belgium, 
Monthly Report on the International Sunspot Number, online catalogue of the 



sunspot index: http://www.sidc.be/sunspot-data/, 2005-2011 



Auger scaler data online, http://auger.colostate.edu/ED/scaler.php 



M. Storini et al., Rome Neutron Monitor, supported by INAF/UNIRomaTre collaboration, 



ht t p : / / cr . izmir an .rssi.ru/ rome / main . htm 



Bray, R.J., Loughhead, R.E. 1979, Sunspots, (Dover Publications, New York) 



Neutron monitor JUNG, http://www.nmdb. eu/?q=node/12 



Neutron monitor APTY, http://www.nmdb. eu/?q=node/118#Apatity 



Fanchiotti, H., Sciutto, S. J., Garcia Canal, C.A., and Hojvat, H. 2004, Fractals, 4, 405 



Peng, C.K., et al. 1994, Phys. Rev. E, 49, 1685 



Hamilton, J. 1994, Time Series Analysis, (Princeton: Princeton Univ. Press) 



-20 - 

MacDonald, D.K.C. 1962 , Noise and Fluctuations: An Introduction, (New York: Wiley) 

Voss R.F., and Clarke, J. 1975, Nature, 258, 317 

Press, W.H. 1978, Comments Astrophys., 7, 103 

Matthaeus W.H. and Goldstein M.L. 1986, Phys. Rev. Lett., 57, 495 

Van Vliet, CM. 1991, Solid-State Electron., 34, 1 

Baillie, R.T. 1996, J. Economet., 73, 5 

Novikov, E., Novikov, A., Shannahoff-Khalsa, D., Schwartz B., and Wright, J. 1997, Phys. 
Rev. E, 56, R2387 

Jensen, H.J. 1998, Self-Organized Criticality, (Cambridge University Press) 

Kantelhardt, J.W., Zschiegner, S.A., Koscielny-Bunde, E., Havlin, S., Bunde, A., and 
Stanley, H.E. 2002, Physica A, 316, 87 

Buldyrev, S.V. 1995, et al, Phys. Rev. E, 51, 5084 

Heneghan, C, and McDarby, G. 2000, Phys. Rev. E, 62, 6103 



This manuscript was prepared with the A AS IATjrpC macros v5.2. 



